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ABSTRACT 


Results  of  experimental  investigations  to  determine  threshold  levels 
for  doubled  neodymium  and  ruby  lasers  are  presented.  The  experimental 
animals  were  rhesus  monkeys. 

Threshold  levels  for  retinal  damage  caused  by  the  doubled  long- 
pulsed  and  Q-switched  neodymium  laser  are  reported,  and  results  of  ex¬ 
periments  to  determine  retinal  damage  and  lens  damage  caused  by  the  Q- 
switched  ruby  laser  are  presented. 

In  addition,  theoretical  considerations  and  calculations  for  retinal 
damage  are  discussed  and  compared  with  experimental  results. 
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FOREWORD 


The  research  program  discussed  in  this  Final  Report  was  conducted 
under  Contract  F41609-70-C-0002,  Project  6301,  Task  630105,  supported 
by  the  USAF  School  oi  Aerospace  Medicine,  Aerospace  Medical  Division, 
Brooks  Air  Force  Base,  Texas,  The  project  monitor  was  MaJ.  I.  Dunsky, 

The  research  reported  covers  the  period  1  September  1069  to  31 
December  1970,  The  project  leader,  Dr.  A,  Vassiliadis,  was  responsible 
for  the  research  activity  under  SRI  Project  8209. 

The  experiments  reported  herein  were  conducted  according  to  the 
"Principles  of  Laboratory  Animal  Care,"  established  by  the  National 
Society  of  Medical  Research. 

Publication  of  this  report  does  not  constitute  Air  Force  approval 
of  the  report's  findings  or  conclusions.  It  is  published  only  lor-  the 
exchange  and  stimulation  of  ideas. 
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ADDENDUM 


The  following  Information  Is  presented  for  the  convenience  of  the  reader, 
to  clarify  certain  aspects  of  the  Final  Report  on  SRI  Project  8209, 
"Ocular  Laser  Threshold  Investigat ions,"  dated  January  1971: 

(1)  The  drugs  used  to  tranquilize  and  anesthetize  the  rhesus 
monkeys  were  the  following: 

•  Sernylan  -  Phencyclidine  Hydrochloride  - 
Parke  Davis  Co.,  Detroit,  Michigan. 

Dose  -  0.04  cc/kg. 

•  Dlabutal  -  Sodium  Pentobarbital  - 
Diamond  Labs.,  Des  Moines,  Iowa. 

Dose  -  0.2  cc/kg. 

(2)  The  focal  lengths  of  lenses  used  to  reduce  the  beam  diameter 
In  the  doubled  Q-switched  neodymium  laser  were  approximately 
6  cm  for  the  positive  lens  and  3  cm  for  the  negative  lens. 

(3)  The  focal  lengths  of  lenses  used  in  the  doubled  long-pulsed 
neodymium  laser  were  approximately  6  cm  for  the  positive  lens 
and  1.7  cm  for  the  negative  lens. 

(4)  The  value  for  the  lens  used  to  correct  for  the  refractive 
error  at  doubled  ruby  wavelength  was  calculated  from 
refractive-index  variation,  as  a  function  of  wavelength, 
estimated  for  the  ocular  media. 

(5)  When  a  monkey  had  an  eye  with  excessive  astigmatism  it  was 
not  used  in  retinal  experiments.  Thus,  all  monkeys  used  in 
retinal  experiments  had  exposures  in  both  eyes. 

(6)  The  titles  of  References  8,  9,  and  10  were  omitted  in  the 
Final  Report  and  are  provided  in  the  following: 

8.  W.  T.  Ham  et  al.,  "Effects  of  Laser  Radiation  on  the 
Mammalian  Eye,"  Trans.  N.Y.  Acad.  Sci.,  Vol.  28,  p.  517 
(1966). 

9.  A.  M.  Clarke  et  al.,  "An  Equilibrium  Thermal  Model  for 
Retinal  Injury  from  Optical  Sources,"  Appl .  Optics, 

Vol.  8,  p.  1051  (1969). 

10.  W.  P.  Hansen  et  al.,  "A  Worst-Case  Analysis  of  Continuous 
Wave  He-Ne  Laser  Hazards  to  the  Eye,"  Appl.  Optics, 

Vol.  6,  p.  1973  (1967). 


I  INTRODUCTION 


Over  the  past  few  years,  concern  for  the  hazards  associated  with 
exposure  to  bright  sources  of  optical  radiation  has  resulted  in  extensive 
studies  by  a  number  of  investigators.  In  this  work,  particular  emphasis 
has  been  given  to  the  damage  effects  of  a  number  of  lasers  in  the  visible 
and  infrared  parts  of  the  spectrum. 

This  report  summarizes  damage  data  for  some  additional  wavelengths 

and  represents  a  continuation  of  work  that  has  been  carried  out  at 

Stanford  Research  Institute  over  the  past  five  years.  Although  data  for 

1  ** 

a  number  of  the  important  lasers  had  been  obtained,  3  there  were  no 
studies  to  cover  the  doubled  output  of  the  two  most  commonly  used  lasers — 
i.e.,  ruby  and  neodymium  doped  materials.  For  this  reason,  the  main 
emphasis  of  the  research  discussed  in  this  report  was  on  damage  thresholds 
for  doubled  neodymium  (530  nm)  and  doubled  ruby  (347  nm)  wavelengths. 

Threshold  experiments  for  doubled  neodymium  were  made  for  both  long- 
pulsed  and  Q- switched  modes  of  operation.  Since  the  doubled  wavelength 
is  transmitted  well  by  the  ocular  media,  retinal  damage  is  of  major  con¬ 
cern,  and  the  threshold  data  were  obtained  for  retinal  damage. 

In  contrast,  the  doubled  output  of  ruby,  being  in  the  ultraviolet, 
is  absorbed  by  the  ocular  media.  Thus,  damage  to  the  cornea,  or  lens  be¬ 
comes  possible.  Experiments  were  therefore  carried  out  to  study  possible 
damage  to  the  lens  as  well  as  to  the  retina. 


References  are  listed  at  the  end  of  the  report. 
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This  report  also  Includes  a  summary  of  some  theoretical  considera¬ 
tions  that  emerged  in  the  course  of  the  program.  The  principal  calcula¬ 
tions  that  are  presented  show,  using  a  simple  model  for  the  retina  and 
the  choroid,  how  the  temperature  variations  oi’  laser  exposures  develop. 
The  theoretical  calculations  are  compared  with  experimental  data,  and 
inconsistencies  are  noted  and  discussed.  It  is  also  argued  that  rate 
processes  in  the  damage  mechanism  must  be  considered  in  order  to  explain 
the  temperature  variations  observed. 

There  is  a  brief  discussion  of  the  rate  processes  in  thermal  damage 
phenomena,  and  an  approximate  estimate  of  the  inactivation  energy  of  the 
proteins  or  enzymes  involved  in  the  threshold-damage  type  of  interaction. 
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II  EXPERIMENTAL  DAMAGE  THRESHOLDS 


A.  General  Remarks 

The  exposures  to  the  retina  and  to  the  lens  were  made  at  a  special 
test  station.  This  delivery  system  is  built  around  a  Zeiss  fundus  camera 
and  has  been  described  in  previous  reports. 3>3  A  photograph  of  the  de¬ 
livery  system  in  use  is  shown  in  Figure  1(a)  and  a  diagram  of  the  system 
is  shown  in  Figure  1(b). 

The  delivery  system  provides  a  means  of  selecting  a  small  part  of 
the  laser  beam  and  accurately  monitoring  it  and  directing  it  into  the 
eye.  In  order  to  select  a  relatively  uniform  part  of  the  laser  beam  an 
adjustable  iris  was  used,  and  the  laser  itself  was  located  at  the  largest 
practical  distance  from  the  delivery  point. 

In  each  of  the  experiments  the  doubled  laser  was  carefully  aligned 
with  the  delivery  system  with  the  aid  of  a  He-Ne  laser.  Calibration  of 
the  delivery  system  was  made  by  placing  a  ballistic  thermopile  at  the 
corneal  plane  of  the  animal's  eye  and  measuring  the  energy  at  the  thermo¬ 
pile  and  at  the  detector  that  was  used  for  monitoring.  This  calibration 
was  made  near  the  laser  so  that  the  whole  doubled  output  could  be  intro¬ 
duced  into  the  thermopile,  thus  providing  adequate  energy  to  perform  the 
calibration. 

The  rhesus  monkeys  used  in  the  experiments  were  given  atropine  for 
at  least  three  days  prior  to  the  day  of  the  exposure.  This  provided 
adequate  paralysis  of  accommodation  in  the  rhesus  eye  so  that  the  experi¬ 
ment  could  be  performed  with  minimum  changes  in  refractive  error.  In 
addition,  this  provided  maximum  dilation  of  the  iris,  which  was  required 
both  for  the  retinal  and  the  lens  exposures. 
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FIGURE  1(a)  PHOTOGRAPH  OF  THE  DELIVERY  SYSTEM 


NARROW 


FIGURE  Kb)  DIAGRAM  OF  THE  DELIVERY  SYSTEM 


A 


The  day  of  the  exposure,  the  bnimals  were  trunquilized  with  Sernylan 
and  then  put  In  deep  anesthesia  with  sodium  pentobarbital.1  At  least  ten 
minutes  were  allowed  to  elapse  from  the  time  of  deep  anesthesia  to  the 
time  the  animal  was  refracted  for  the  experiment.  This  assured  that  the 
changes  in  refraction  that  follow  the  drug  administration  had  subsided. 
The  refraction  was  made  to  the  nearest  0.25  diopter  and  the  appropriate 
lens  was  placed  in  front  of  the  animal's  eye  during  the  experiment.  If 
an  animal  had  astigmatism  of  over  0.5  diopter  it  was  not  used  in  the 
retinal  experiments. 

B.  Doubled  Q-Switched  Neodymium 
1.  Laser  Characteristics 

The  neodymium  laser  cavity  comprised  a  6- inch  rod  1/2  inch  in 
diameter  that  has  Brewster-angle  cuts.  The  rod  was  water-cooled  an'i 
pumped  with  two  air-cooled  linear  flash  lamps.  The  laser  was  Q-switched 
by  a  spinning  quartz  prism  turning  at  500  rps  that  was  optimumly  timed 
to  switch  the  cavity  Q.  A  sapphire  flat  Fabry-Perot  served  as  the  output 
end. 


The  output  beam  was  reduced  in  diameter  by  a  factor  of  two  by 
the  use  of  two  lenses.  The  beam  was  then  passed  through  an  accurately 
oriented  ADP  crystal  with  which  the  530-nm  radiation  was  obtained.  A 
photograph  of  a  typical  pulse  is  shown  in  Figure  2.  The  conversion 
efficiency  in  the  doubler  was  not  very  good— typically  3  to  5  percent. 
However,  this  provided  more  than  adequate  energies  at  the  test  station, 
as  described  in  the  next  section. 

The  doubled  Q-switched  laser  was  located  approximately  35  feet 
from  the  delivery  system.  The  divergence  of  the  beam  at  the  test  station 
was  measured  to  be  1.7  mrad.  Thus,  the  spot  size  at  the  retina,  although 
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FIGURE  2  TYPICAL  PULSE  OF  DOUBLED  O-SWITCHED  NEODYMIUM  LASER 


determined  by  the  animal's  eye  optics,  w 1 1 1  be  greater  than  approximately 
23  ^ 


2.  Experiments  and  Results 

The  experiimnts  initially  involved  a  large  number  ol  exposures 
in  the  paramacular  area.  In  initial  expert  merits,  animals  received  16 
paramacular  exposures  and  lour  macular  exposures.  The  paramacular  ex¬ 
posures  were  all  placed  within  the  temporal  superior  and  interior  vessel 
arcades.  Thus,  the  exposures  were  m  the*  immediate  neighborhood  ol  the 
macula  and  near  the  posterior  pole.  Since  the  primary  interest  is  in  the 
macular  threshold  level,  the  large  majority  ol  the  experiments  were  made 
with  eight  paramacular  exposures  and  tour  macular  exposures  in  each  eye. 

The  procedure  used  was  similar  to  the  protocols  ol  previous 
experiments,*"’^  Thus,  the  laser  exposures  were  placet!  in  the  eye  and 
the  observation  to  determine  whether  a  lesion  has  lortued  was  made  one 
hour  alter  exposure.  Each  exposure*  site  was  examined  and  classilied  as 
either  damaged  or  not  damaged. 
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The  retinal  response  for  these  exposures  is  similar  to  those 
observed  in  previous  Q-swltched  experiments,1*3  Thus,  near  threshold, 
no  immediate  response  is  seen.  As  the  level  is  increased,  however,  a 
point  is  reached  where  an  immediate  response  is  seen.  This  has  the 
appearance  of  a  small  white  patch  interpreted  as  edema,  easily  visible 
through  an  ophthalmoscope. 

The  experimental  results  presented  in  this  section  are  a  summary 
of  the  data  on  20  rhesus  monkeys.  The  data  are  presented  separately  as 
macular  and  paramacular.  For  each  set  of  data  the  individual  exposures 
are  shown  in  a  figure,  and  a  probability  curve  is  presented  for  each 
case. 

The  data  for  the  paramacular  area  are  summarized  in  Figure  3. 
From  the  data  it  is  apparent  that  the  low  side  of  the  damage  level  was 
stressed  in  order  to  assure  that  the  threshold  that  is  obtained  is 
meaningful.  The  data  show  that  the  50-percent  probability  for  damage 
occurs  at  about  10  p.J.  This  may  be  compared  with  the  Q-switched  ruby 
data,  which  show  a  50-percent  damage  probability  at  about  20  ^J.1 

The  macular  data  summarized  in  Figure  4  show  a  50-percent 
probability  for  damage  at  about  5  m.J.  This  is  appreciably  lower  than 
for  the  paramacular  area.  This  same  trend  has  been  observed  for  Q- 
switched  neodymium,3  long-pulsed  ruby,1  and  short  He-Ne  laser  exposures.4 

C.  Doubled  Long-Pulsed  Neodymium 
1.  Laser  Characteristics 

The  neodymium  laser  cavity  used  in  this  experiment  was  the  same 
as  that  used  in  the  Q-switched  experiments.  Instead  of  a  spinning  prism, 
however,  a  stationary  homosi 1-quartz  roof  prism  was  used  as  the  back  end 
of  the  cavity.  The  laser  was  upo rated  in  a  long-pulse  mode  with  a  typical 
pulse  train  lasting  approximately  600  ^s,  as  shown  in  Figure  5. 
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FIGURE  3  PARAMACULAR  DAMAGE  TO  RHESUS  MONKEY  RETINA  WITH  DOUBLED  Q-SWITCHED  NEODYMIUM 


FIGURE  4  MACULAR  DAMAGE  TO  RHESUS  MONKEY  RETINA  WITH  DOUBLED  Q-SWITCHED  NEODYMIUM 


TA-8209-71 


FIGURE  5  TYPICAL  PULSE  TRAIN  OF  LONG-PULSED  NEODYMIUM  LASER 

The  output  bean  was  reduced  in  diameter  by  a  factor  of  about 
3.5  by  the  use  of  two  lenses,  in  order  to  obtain  higher  power  density 
and  thus  increase  conversion  efficiency.  The  beam  was  then  passed 
through  an  ADP  crystal,  with  which  the  530  nm  radiation  was  obtained. 

The  conversion  efficiency  in  the  doubler  was  not  very  good.  Considerable 
difficulty  was  experienced  in  optimizing  the  doubled  output.  It  was 
found  necessary  to  incorporate  an  accurate  thermal  control  over  the  laser 
rod.  With  a  system  that  maintained  the  rod  to  within  0.6°  C  it  was 
possible  to  obtain  reproducible  results,  and  the  crystal  could  be  tuned. 

A  compromise  was  established  between  the  point  of  maximum  output  and  the 
point  where  variations  of  temperature  had  minimal  effect  on  the  output. 

The  doubled  laser  was  located  approximately  12  feet  away  from 
the  delivery  system.  This  was  necessary  in  order  to  obtain  sufficient 
energy  in  the  eye  to  cause  minimal  damage. 

Even  though  the  laser/doubler  combination  was  close  to  the 
delivery  station,  the  beam  was  relatively  well  collimated,  because  of 
the  critical  phase-matching  requirement  in  the  doubling  crystal.  The 
beam  divergence  was  measured  to  be  2  mrad.  Thus,  the  spot  size  in  the 
retina  for  rhesus  monkeys  would  be  larger  than  27  depending  on  the 
eye  optics. 
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2.  Experiments  and  Hesults 


Because  of  the  difficulties  In  obtaining  sufficient  doubled 
output  and  the  fact  that  the  doubling  crystal  wus  being  damaged,  it  was 
decided  to  keep  the  number  of  exposures  to  a  minimum.  The  paramacular 
exposures  were  all  placed  within  the  temporal  superior  and  inferior 
vessel  arcades.  Thus,  the  exposures  were  in  the  immediate  neighborhood 
of  the  macula  and  near  the  posterior  pole.  Since  the  primary  interest 
Is  In  the  macular  threshold  level,  there  were  four  paramacular  exposures 
and  four  macular  exposures  in  each  eye. 

The  procedure  used  was  similar  to  that  used  in  the  Q-switched 
case  discussed  above.  The  appearance  of  the  retinal  lesions  at  threshold 
and  above  threshold  is  similar  to  that  reported  previously  for  other 
lasers.3  >3 

The  experimental  results  presented  in  this  section  are  a  summary 
of  the  data  on  20  rhesus  monkeys.  The  data  are  presented  separately  as 
macular  and  paramacular.  The  data  for  the  paramacular  area  are  summarized 
in  Figure  6.  The  data  show  that  the  50-percent  probability  for  damage 
occurs  at  about  74  jiJ.  The  macular  data  summarized  in  Figure  7  show  a 
50-percent  probability  for  damage  of  about  38  pj.  This  is  appreciably 
lower  than  the  paramacular  area.  This  same  trend  has  been  observed  for 
a  number  of  the  lasers3-4  as  well  as  for  doubled  Q-switched  neodymium 
discussed  in  the  previous  section. 

D.  Doubled  Q-Switched  Ruby 

1.  Laser  Characteristics 

The  laser  makes  use  of  a  6- inch- by- 1/2- inch  ruby  rod  with 
Brews ter-angle  cuts.  The  rod  is  pumped  by  two  12-mm-bore  linear  flash 
lamps  in  a  close-coupled  geometry.  The  lamps  are  air-cooled,  whereas 
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FIGURE  7  MACULAR  OAMAGE  TO  RHESUS  MONKEY  RETINA  WITH  DOUBLED  LONG-PULSED  NEODYMIUM 


the  ruby  rod  is  water- coo  led.  The  laser  is  Q- switched  with  a  spinning 
quartz  prism  revolving  at  300  rps.  Typically,  the  out  put -pulse  energy 
of  the  ruby  laser  at  60  1  nm  is  1  .]  and  tlie  pulse  length  is  30  ns 

The  output  ol  the  laser  \^as  introduced  immediately  into  a 
properly  oriented  ADP  crystal  winch  doubled  the  ruby  ircquency,  and  a 
pulse  was  obtained  at  317  nm.  Immediately  alter  the  doubler  a  lilter 
was  used  to  pass  only  the  doubled  output  and  absorb  the  lundamental  ruby 
wavelength.  A  photograph  ol  a  typical  pulse  is  shown  m  figure  H. 
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FIGURE  8  TYPICAL  OUTPUT  PULSE  OF  Q-SWITCHED  RUBY  LASER 


The  doubled  laser  was  located  33  ieet  1  rom  the  delivery  system. 
The  beam  divergence  measured  lor  this  laser  was  1.5  mrad.  Thus,  when 
properly  focused  in  the  rhesus  monkey  eye  the  retinal  spot  size  would  be 
larger  than  20  fA,  depending  on  the  eye  optics. 

For  the  lens-damage  experiments,  the  size  ol  the  loeused  spot 
provided  by  the  quartz  lens  that  was  used  was  estimated  to  be  HO  _  in 
diameter. 

2.  Retinal  Experiments  and  Results 

The  doubled  ruby  laser  was  used  in  a  retinal  threshold  study. 
The  exposures  in  rhesus  eyes  were  all  placed  within  the  temporal  superior 
and  inferior  arcades.  Thus  all  exposures  were  either  in  the  macula  (>r 
in  the  immediate  neighborhood  of  the  macula. 


The  procedure  was  similar  to  that  used  in  previous  experiments. 
Thus,  the  criterion  for  damage  was  the  requirement  for  an  ophthalmoscopi- 
cally  visible  response  within  one  hour  after  exposure. 

A  total  of  eight  rhesus  monkeys  were  used  for  these  exper iments-- 
i.e.,  a  total  of  sixteen  eyes.  In  half  of  the  eyes,  the  exposures  were 
made  with  the  proper  correction  lens  located  in  front  of  the  eye.  In 
the  other  half  of  the  eyes,  a  negative  lens  of  approximately  2.5  diopters 
was  included  in  front  of  the  eye.  This  lens  was  used  to  correct  for  the 
refractive  power  of  the  eye  at  347  nm,  which  is  much  higher  than  the 
value  obtained  with  retinoscopy,  as  estimated  from  calculations. 

In  all  these  exposures,  no  retinal  damage  was  observed.  The 
data  are  summarized  in  Figure  9,  where  it  is  seen  that  up  to  the  highest 
energy  values  that  were  used — i.e.,  somewhat  over  6  mJ--no  visible  re¬ 
sponse  was  observed  in  the  retina. 

3.  Lens-Damage  Experiments  and  Results 

The  doubled  ruby  laser  was  also  used  to  study  lens  damage  in 
the  rhesus  monkey  in  order  to  determine  a  threshold  level. 

The  rhesus  monkey  eyes  were  dilated  maximally  by  application 
of  atropine  prior  to  the  day  of  the  experiment.  The  laser  beam  was  then 
focused  by  the  use  of  a  quartz  lens  of  approximately  2- inch  focal  length. 
The  focused  beam  was  aimed  at  the  edge  of  the  lens  near  the  iris  at  four 
well  separated  locations.  Thus,  each  eye  received  four  exposures  in  the 
lens.  Each  eye  was  examined  carefully  under  a  slit-lamp  immediately 
after  exposure  and  at  one  month  and  at  three  months  post-exposure. 

A  total  of  25  animals  were  used  in  this  study  of  lens  damage, 
and  the  results  discussed  below  are  a  summary  of  these  exposures. 
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FIGURE  9  SUMMARY  OF  RETINAL  EXPOSURES  OF  RHESUS  MONKEYS  TO  DOUBLED 
Q-SWITCHED  RUBY 


The  lenticular  damage  observed  consisted  of  small  areas  of 
opacification — initially  very  tenuous  and  difficult  to  photograph.  At 
the  end  of  one  month,  however,  the  opacities  appeared  considerably  more 
well  defined  and  in  some  cases  considerably  larger  than  was  thought  to 
be  the  exposed  area.  The  appearance  of  the  opacities  in  detail  show 
small  white  granules  dispersed  throughout  the  damaged  area.  Some  ex¬ 
amples  of  the  opacities  taken  at  one  and  three  months  post-exposure  are 
shown  in  Figures  10  through  13.  It  will  be  noted  that  there  is  con¬ 
siderable  difference  in  the  sizes  of  the  lesions.  It  must  be  noted  that 
all  these  are  considerably  above  the  threshold  level. 

The  data  based  on  the  immediately  post-operative  observation 
are  summarized  in  Figure  14.  The  data  for  3-month  post-exposure  examina¬ 
tion  are  summarized  in  Figure  15.  Statistically,  the  threshold  level  did 


16 


FIGURE  10  SMALL  LENS  OPACITIES  IN  RHESUS  MONKEY  ONE  MONTH  AFTER 
EXPOSURE  TO  DOUBLED  Q-SWITCHED  RUBY 
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FIGURE  11 


RHESUS  MONKEY  LENS  OPACITIES  ONE  MONTH  AFTER  EXPOSURE 
TO  DOUBLED  Q-SWITCHED  RUBY 
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FIGURE  12 


RHESUS  MONKEY  LENS  OPACITY  THREE  MONTHS  AFTER  EXPOSURE 
TO  DOUBLED  Q-SWITCHED  RUBY 
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FIGURE  13 


POSTERIOR  SURFACE  LENS  CHANGE  IN  RHESUS  MONKEY  AFTER 
DOUBLED  RUBY  LASER  EXPOSURE 
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FIGURE  14  RHESUS  MONKEY  LENTICULAR  DAMAGE  BASED  ON  IMMEDIATE  POST-OPERATIVE  OBSERVATION 
FOR  EXPOSURE  TO  DOUBLED  Q-SWITCHED  RUBY 
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not  change  appreciably.  Thus,  although  some  damage  sites  that  were  seen 
nt  Inter  examination  were  not  seen  immediately,  there  were  also  some  that 
were  not  seen  in  later  examinations  although  they  were  marginally  observed 
earl  ier . 

It  is  noted  that  the  50-percent  probability  for  damage  in  these 

2 

experiments  was  approximately  0,7  mJ.  This  is  equivalent  to  14  J/cm  at 
the  site  of  the  exposure. 
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Ill  THEORETICAL  CONSIDERATIONS  OF  RETINAL  DAMAGE 


A.  Absorption  of  Light  in  the  Fundus 

In  the  visible  and  near  infrared  parts  of  the  spectrum  the  ocular 
media  are  relatively  transparent,  and  thus  light  is  transmitted  effec¬ 
tively  to  the  back,  or  fundus,  of  the  eye.  The  light  arriving  at  the 
retina  at  a  high  level  of  illumination  passes  through  the  various  layers 
of  the  retina  virtually  unattenuated  so  that  a  high  percentage  of  the 
light  incident  on  the  cornea  reaches  the  pigment  epithelium,  the  first 
absorbing  layer  in  the  fundus.  The  highest  absorption  per  unit  volume 
takes  place  in  this  layer. 

The  high  absorption  of  the  pigment  epithelium  layer  is  due  to  the 
high  concentration  of  melanin  granules  that  are  located  in  the  pigment 
epithelium  cells.  These  granules  are  ellipsoids  and  spheres  in  the  range 
of  1/2  to  1  p  in  diameter.  It  is  usually  noted  that  these  granules  are 
located  in  the  anterior  part  of  the  pigment  epithelium  cells.  Thus, 
although  the  pigment  epithelium  cells  are  approximately  10  p,  in  thick¬ 
ness,  the  melanin  granules  are  located  within  a  3-to-4-p-thick  band 
immediately  adjacent  to  the  outer  extremities  of  the  rods  and  cones. 

An  appreciable  fraction  of  the  light  is  absorbed  by  the  pigment 
epithelium.  The  majority,  however,  proceeds  to  the  choroid,  where  it  is 
scattered  and  partly  absorbed  by  the  pigment  granules  and  hemoglobin  are 
present  in  the  choroid.  The  light  that  is  not  absorbed  by  the  choroid 
proceeds  to  the  sclera. 
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B.  Thermal  Damage 

Although  the  processes  associated  with  vision  are  photochemical  In 
nature,  the  retinal  damage  caused  by  intense  sources  of  light  Is  generally 
understood  to  be  thermal  In  nature.  This,  at  least,  appears  to  be  true 
for  damage  associated  with  acute  exposures  lasting  up  to  tens  of  seconds. 
Thus,  the  light  that  Is  absorbed  by  the  pigment  epithelium  and  choroid 
Is  converted  Into  heat.  Consequently,  If  the  intensity  is  high  enough 
the  temperature  rise  of  the  absorbing  tissues  can  be  sufficient  to  cause 
thermal  damage. 

The  normal  operation  of  a  living  cell  requires  that  the  ambient 
temperature  be  kept  within  very  restrictive  limits.  It  is  not  surprising, 
therefore,  that  a  relatively  small  increase  in  temperature  can  lead  to 
damaging  consequences.  Many  theories  have  been  proposed  as  to  how  thermal 
damage  of  the  cell  takes  place.  Among  the  current  theories  are  that 
damage  is  due  to  alterations  of  nucleic  acids,  inactivation  of  enzymes, 
or  denaturation  of  proteins.5 

These  inactivation  and  denaturation  processes  must  be  analyzed  as 
rate  processes.  Just  as  heat  speeds  up  chemical  reactions,  a  small 
temperature  rise  can  profoundly  Influence  the  rate  of  the  inactivation 
or  denaturation  processes.  Thus,  In  model  calculations,  some  considera¬ 
tion  must  be  given  not  only  to  the  thermal  history  of  the  tissues  but 
also  the  rate  processes  that  are  involved. 

C.  Models  for  Retinal  Damage 

The  first  attempt  at  a  theory  for  retinal  burns  was  developed  by 
Vos.5  He  assumed  a  thin  (10-^- thick)  rectangular  region  absorbing  uni¬ 
formly,  and  calculated  the  thermal  history  of  the  exposed  site.  He  con¬ 
sidered  times  that  were  longer  than  20  ms  and  spot  sizes  that  were  greater 
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than  180  jx.  Vob  concluded  thot  the  retinal  damage  Indicated  by  data  of 
Ham  et  al.7  was  closely  associated  with  the  calculated  amount  of  steam 
formation  produced  by  the  exposure.  Hc'did  point  out,  however,  that 
damage  should  occur  at  lower  levels  of  irradiation  since  albumins  begin 
to  coagulate  at  temperatures  between  45°  C  and  60“  C. 

Subsequently,  Ham  ct  al.s  developed  a  model  in  which  heat  was 
generated  uniformly  in  the  pigment  epithelium  and  in  the  choroid.  More 
recently  Clarke  et  al.*  reported  on  a  model  in  which  pigment  epithelium 
and  choroid  were  assumed  to  act  as  uniform  absorbers.  In  the  latter 
case  only  a  steady-state  model  was  treated.  Another  analysis  of  thermal 
damage,  also  for  steady-state  condition,  was  made  by  Hansen  et  al.,10 
who  considered  a  steady-state  solution  of  a  cylindrical  volume  of  uniform 
heat  generation  representing  the  pigment  epithelium. 

All  these  models  have  shortcomings  and  difficulties,  and  comparison 
with  experimental  data  has  not  been  satisfactory.  None  of  the  models 
takes  into  account  a  rate  equation  for  the  denaturation  or  inactivation 
process.  In  the  following  sections  we  consider  a  simple  model  in  order 
to  make  calculations  to  compare  with  recent  experimental  damage  data  for 
CW  lasers. 

D.  Temperature  Variation  in  Retinal  Exposures 

In  this  section  we  consider  a  simple  model  for  calculating  the 
temperature  rise  at  the  center  of  a  retinal  exposure  to  intense  radiation. 
We  consider  the  time  regime  in  which  thermal  relaxation  plays  a  strong 
role.  These  arguments  would  thus  be  most  appropriate  for  exposures  to 
CW  lasers. 

We  assume  that  the  beam  that  is  incident  on  the  retina  has  variations 
only  in  the  radial  direction.  Two  cases  are  considered:  first,  the 
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intensity  is  assumed  uniform  und  of  radius  u,  and  second,  a  gaussian  dis¬ 
tribution  is  assumed. 

The  absorption  is  assumed  to  take  place  in  two  layers,  the  first 
representing  the  pigment  epithelium  and  the  second  representing  the 
choroid.  It  is  of  interest  then  to  calculate  the  temperature  on  the 
axis,  since  the  highest  temperature  will  occur  somewhere  on  the  axis. 

For  any  point  on  the  axis  the  temperature  rise  con  be  calculated  as  the 
sum  of  the  contributions  from  the  pigment  epithelium  layer  and  from  the 
choroid. 


As  shown  in  the  Appendix,  the  solution  for  an  absorbing  layer  of 
thickness  2 £  is  given  by  Eq.  (A-12)  as 


for  a  uniform  beam  of  radius  a.  For  a  beam  of  gaussian  distribution  it 
is  shown  by  Eq.  (A-14)  that  the  temperature  rise  is  given  by: 
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where  the  radial  distribution  is  given  by 
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In  Eqs.  (1)  and  (2), 


and 


k  =  Diffusivity  (cm  /s) 
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whoro 


=  Power  density  (wntts/cm  ) 

2 

K  =  Heat  conductivity  (cal  cm/cm  -s-  C) 


0  =  Fraction  of  power  ubsorbcd. 
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computer  calculations  were  made  for  | , 
dimensionless  variables  dcf  ined  by 
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Results  were  obtained  for  a  number  of  parameter  variations  so  that 
various  spot  sizes  and  choroid  and  pigment  epithelium  configurations  could 
be  selected.  Some  of  the  results  are  plotted  in  Figures  16  and  17  for 
the  uniform  beam  cross  section,  and  in  Figures  18  and  19  for  the  gaussian 
beam  cross  section. 


Before  applying  the  simple  thermal  model  to  experimental  configura¬ 
tions  so  that  a  comparison  with  the  empirical  data  can  be  made,  thermal 
damage  of  biological  material  is  discussed  in  the  next  section. 


E.  Thermal  Rate  Processes 

The  denaturation  of  proteins  and  the  inactivation  of  enzymes  must 
be  analyzed  as  a  rate  process.  As  shown  by  observations  of  rates  of 
chemical  reactions,  the  empirical  expression  for  the  reaction  rate  con¬ 
stant  has  the  form 
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FIGURE  16  TEMPERATURE  RISE  AS  A  FUNCTION  OF  TIME  FOR  VARIOUS  SPOT  SIZES  AT  THE  CENTER  OF 
AN  ABSORBING  DISC  —  UNIFORM  BEAM  CROSS  SECTION 
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FIGURE  17  TEMPERATURE  RISE  AS  A  FUNCTION  OF  TIME  FOR  VARIOUS  DISTANCES  FROM  CENTER  OF  A  DISC 
UNIFORM  BEAM  CROSS  SECTION 
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DISC  —  GAUSSIAN  BEAM  CROSS  SECTION 
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OISC  —  GAUSSIAN  BEAM  CROSS  SECTION 
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FIGURE  19  TEMPERATURE  RISE  AS  A  FUNCTION  OF  TIME  FOR  VARIOUS  DISTANCES  FROM  CENTER  OF  DISC 
GAUSSIAN  BEAM  CROSS  SECTION 
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This  can  also  be  written  ns 


k  =  A  e 
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In  a  similar  way,  from  statistical  mechanics  treatments11  it  can  be 

shown  that  the  reaction  rote  is  given  by 

*  *  * 
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where  k  is  the  Boltzmann  constant,  h  is  Planck's  constant,  AF  is  the 

* 

change  in  activation  free  energy,  AH  the  change  in  heat  energy,  and 
* 

AS  the  change  in  entropy. 


This  may  be  rewritten  as 
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AH  +  RT 
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RT 
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which  is  seen  to  be  similar  to  the  Arrhenius  equation,  Eq.  (6),  if  we 
identify  with  (AH*  +  RT) .  Since  the  influence  of  temperature  on  the 
reaction  rate  is  dominated  by  the  exponential  factor,  although  A  (Eq.  7) 
does  vary  slightly  with  temperature,  it  is  usual  to  assume  that  A  is 
independent  of  temperature  in  the  treatment  of  experimental  data. 

The  reaction  associated  with  denaturation  of  proteins  and  inactiva¬ 
tion  of  enzymes  is  usually  unimolecular 2  thus  the  reaction  follows 
first-order  kinetics,  and  we  can  write 
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where  C  is  the  concentration. 

The  surviving  traction  S  may  then  be  expressed  as 


(11) 


(12) 


In  terms  of  the  temperature  rise  in  tissues  t  we  therefore  must 
consider  the  equation 
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where  T  is  the  ambient  temperature  and  4r  is  the.  temperature  rise  due  to 
o 

the  radiation  being  absorbed  by  the  medium. 


In  comparing  the  theoretical  results  with  the  experimental  data  it 

must  be  kept  in  mind  that  the  experimental  results  usually  involve  a 

gross  examination  of  the  exposed  site  and  it  is  unknown  what  the  survival 

fraction  is  when  damage  is  observed.  However,  since  the  solution  of  Eq. 

(13)  is  relatively  insensitive  to  the  survival  fraction  that  is  used,  we 

can  select  arbitrarily  a  value  for  C/C  =  1/e. 

0 

Thus,  the  solution  of  Eq.  (13)  in  the  form 


RfT-HjKt)] 

A  e  dt  =  1  (14) 
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must  bo  obtained  for  the  temper  a  lure  rise,  -Jj ,  given  by  the  thermal  con¬ 
duction  model. 

F .  Comparison  of  Theory  and  Experimental  Results 

in  this  section  wo  consider  the  application  of  simple  thermal  models 
to  the  understanding  of  experimental  retinal  threshold  data  obtained  with 
CW  lasers.  The  experimental  threshold  data  that  are  examined  were  ob¬ 
tained  using  an  argon  CW  laser  in  rhesus  monkey  eyes. 

The  data  shown  in  Figure  20  represent  threshold  data  for  an  argon 
laser.3  Each  point  on  the  plot  represents  the  50-percent  probability 
for  damage  value,  and  is  a  result  of  a  large  number  of  exposures  in  rhesus 
monkey  eyes.  The  data  accumulated  were  based  on  an  ophthalmoscopically 
visible  criterion.  Two  features  of  the  data  shown  in  Figure  20  are  of 
interest.  First  it  is  noted  that  the  curves  continue  to  decrease  as  a 
function  of  the  length  of  the  exposure.  This  trend  is  much  more  pro¬ 
nounced  than  would  be  expected  from  the  constant- temperature  thermal- 
damage  model  that  is  usually  referred  to  in  thermal-damage  discussions. 

The  second  feature  noted  in  the  data  shown  in  Figure  20  is  that  the 
spacing  between  the  curves  does  not  follow  the  simple  relationship  ex¬ 
pected  from  simple  model  considerations.  Thus,  the  data  for  the  different 
spot  sizes  fall  much  closer  together  than  would  be  predicted  from  simple- 
model  arguments. 

We  turn  now  to  a  discussion  of  the  first  feature — the  shape  of  the 
damage- threshold  curves.  We  consider  the  data  for  the  30-^  retinal  spot, 
and  present  calculations  for  a  number  of  simple  models  to  indicate  the 
trends  observed. 

First,  a  simple  one-layer  absorbing  model  is  used  to  simulate  the 
pigment  epithelium  (PE)  absorption.  Thus,  the  beam  is  assumed  to  be 
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absorbed  only  in  Ihc  pigment  epithelium,  and  the  temperature  at  the 
center  ol  the  PE  is  calculated.  The  PE  luyer  is  assumed  to  be  4  ^  in 
thickness,  which,  it  is  felt,  realistically  represents  the  absorbing 
part  of  the  layer.  Using  the  curves  of  Figure  16  and  Eqs.  (4)  and  (5) 
we  find  that  the  threshold  variation  for  this  simple  model  must  have  the 
form  shown  in  Figure  21.  The  curve  has  been  arbitrarily  shown  to  match 


TIME  —  ms 


T  A -8209 -84 


FIGURE  21  COMPARISON  OF  RELATIVE  SHAPE  OF  THRESHOLD  VARIATION  FOR  SIMPLE 
MODEL  WITH  ACTUAL  DATA  —  ASSUMING  ONLY  PIGMENT  EPITHELIUM 
LAYER  ABSORBING  AND  PE  AND  CHOROID  EQUALLY  ABSORBING 


the  experimental  data  at  one  second.  It  is  apparent  that  the  experimental 
data  show  considerably  more  change  with  time  than  would  be  expected  from 
simple  theory.  Very  little  variation  is  to  be  noted  between  a  4-^  and  a 
10-|i  thickness  (10  ^  is  the  thickness  usually  assumed  for  the  PE  layer). 

A  comparison  between  the  two  is  shown  in  Figure  22. 
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FIGURE  22  COMPARISON  OF  RELATIVE  SHAPE  OF  THRESHOLD  VARIATION  FOR  TWO  THICKNESSES  OF  THE 
PIGMENT  EPITHELIUM  LAYER  FOR  VARIOUS  RETINAL  SPOT  SIZES 


FIGURE  22  COMPARISON  OF  RELATIVE  SHAPE  OF  THRESHOLD  VARIATION  FOR  TWO  THICKNESSES  OF  THE 
PIGMENT  EPITHELIUM  LAYER  FOR  VARIOUS  RETINAL  SPOT  SIZES 


A  more  complex  model  wuh  then  assumed  so  an  to  represent  the  PE  and 
the  choroid  more  real  1st  Icnlly.  It  is  assumed  that  the  choroid  absorption 
layer  is  *10  thick  and  is  centered  30  ^  behind  the  pigment  epithelium 
absorbing  layer.  It  is  then  assumed  that  an  equal  amount  of  power  is 
absorbed  by  the  PE  and  by  the  choroid,  A  calculation  of  the  threshold 
power  lor  constant  temperature  of  this  model  is  also  shown  in  Figure  21, 

Here  it  is  seen  that  the  curve  is  appreciably  closer  to  the  actual  ex¬ 
perimental  curve;  however,  it  Is  still  insufficient  to  explain  the  ex¬ 
perimental  observations.  It  is  to  be  noted  that  if  the  large  majority 
of  the  power  was  absorbed  by  the  choroid,  the  calculated  curve  could 
come  closer  to  the  experimental  observations.  This,  however,  would  be 
contrary  to  absorption  measurements  of  the  various  absorbing  layers. 

We  turn  now  to  thermal-damage  calculations  that  model  the  damage 
process  by  considering  the  rate  processes  and  by  solving  Eq.  (14)  for 
the  temperature- rise  variation  prescribed  by  the  thermal-model  calculations. 

It  was  assumed  that  the  fundus  consists  of  two  absorbing  layers, 

representing  the  PE  and  choroid,  that  absorb  equal  amounts  of  power. 

Equation  (14)  was  then  used  to  calculate  requirements  of  exposure  time 

for  a  number  of  inactivation  energies,  E^,  as  a  function  of  input-power 

levels.  Figure  23  shows  examples  for  E  =  50,000,  70,000,  and  100,000 

cal/mole,  which  represents  a  range  of  values  typical  of  low  inactivation 

energies  lor  proteins  and  enzymes.  A  visual  fit  of  these  curves  with  the 

experimental  data  showed  that  the  curve  for  E  =  70,000  cal/mole  gave  a 

good  fit.  This  is  shown  in  Figure  24,  where  the  agreement  is  seen  to  be 

excellent.  It  must  be  kept  in  mind,  however,  that  this  match  depends  on 

the  absorbing  model  that  is  chosen — i.e.,  the  relative  absorption  between 

the  PE  and  the  choroid  influences  the  value  of  E  that  is  calculated. 

o 

Thus,  a  more  accurate  knowledge  of  the  absorption  of  the  PE  and  choroid 
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FIGURE  23  RELATIVE  THRESHOLD  VARIATIONS  FOR  THREE  INACTIVATION  ENERGIES 
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FIGURE  24  COMPARISON  OF  RELATIVE  SHAPE  OF  THRESHOLD  VARIATION  FOR  INACTIVATION  THERMAL 

MODEL  WITH  EXPERIMENTAL  DATA.  EQ  =  50.000  cal/mole.  PE  AND  CHOROID  ASSUMED  EQUALLY 
ABSORBING. 


would  bo  necessary  beloro  un  accurate  value  of  E^  could  be  determined 
and  an  identification  of  «  particular  enzyme  or  protein  could  be  made. 

Finally,  we  consider  the  problem  oi  variation  of  threshold  with 
spot  size.  From  the  calculations  of  temperature  shown  in  Figures  16 
through  19  we  obtain  vulucs  ior  the  power  level  required  to  reach  the 
same  temperature  at  each  spot  size,  at  a  time  of  80  ms,  where  we  have 
three  experimental  data  points.  Calculations  were  made  for  the  simple 
PE  model  as  well  as  the  two-layer  PE  and  choroid  model.  The  relative 
power  levels  required  are  shown  in  Figure  25,  where  for  each  curve  the 
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FIGURE  25  COMPARISON  OF  SIMPLE  MODELS  —  RELATIVE 
SHAPE  OF  THRESHOLD  VARIATION  AS  A 
FUNCTION  OF  SPOT  SIZE 
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normalization  i s  made  arbitrarily  to  the  200-^  spot  size.  It  is  apparent 
that  there  is  eons tderable  discrepancy  between  the  theory  and  the  experi¬ 
mental  data.  Although  the  addition  of  the  choroid  does  bring  tho  calcu¬ 
lation  closer  to  the  experiment,  choroid  absorption  much  higher  than  the 
PE  absorption  would  be  required  to  bring  the  calculation  and  the  experi¬ 
mental  data  closer  together. 

Application  of  the  inac t ivat ion- rate-process  calculations  would  not 
move  the  theoretical  curves  appreciably  closer  to  the  experimental  points; 
thus  there  appears  to  be  a  need  for  u  nonthcrmal  explanation  of  this 
phenomenon.  As  mentioned  earlier,3*1 3  the  most  reasonable  explanation 
may  be  that  the  effect  is  due  to  the  scattering  of  the  ocular  refracting 
surfaces  and  media. 
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IV  CONCLUSIONS 


Threshold  levels  for  retinal  damage  were  reported  for  rhesus  monkeys 
for  a  doubled  Q-switched  neodymium  laser.  The  threshold  was  found  to  be 
10  n«J,  which  is  lower  than  for  Q-switched  ruby.  The  macular  area  was 
found  to  be  approximately  twice  as  sensitive  as  the  paramacular  area. 

Threshold  data  for  doubled  long-pulsed  neodymium  were  taken,  and  a 
threshold  of  74  was  found  to  be  approximately  seven  times  higher  than 
for  the  Q-switched  mode. 

Threshold  data  for  retinal  damage  for  the  doubled  Q-switched  ruby 
laser  produced  no  visible  damage  up  to  levels  of  6  mj.  Thus,  the  absorp¬ 
tion  by  the  ocular  media  does  increase  the  threshold  at  347  mfi  considerably. 

2 

Threshold  data  on  lenticular  damage  gave  a  value  of  approximately  14  J/cm 
as  the  required  incident  level  to  cause  damage.  This  is  not  to  say  that 
lens  damage  will  occur  before  retinal  damage  occurs  at  this  wavelength. 

Theoretical  models  for  retinal  damage  were  discussed  and  computer 
parametric  calculations  were  presented.  Comparison  of  calculated 
threshold  values  with  experimental  data  has  shown  some  disagreement  with 
simple  thermal  models.  It  was  shown  that  the  experimental  variation  of 
threshold  with  spot  size  is  inconsistent  with  simple  thermal  models.  In 
addition,  it  was  shown  that  the  experimentally  observed  variation  of 
threshold  with  time  could  not  be  explained  by  simple  thermal  models. 

However,  it  was  shown  that  the  time  variations  observed  in  the  experi¬ 
mental  data  are  consistent  with  a  theory  that  describes  the  denaturation 
or  inactivation  process  of  the  damage  mechanism  with  an  appropriate  rate 
equation. 
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Appendix 


HEAT- FLOW  PROBLEM 

The  conduction  of  heat  in  un  Isotropic  solid  Is  described  by  the 
following  equation: 


& 

at 


2 

hv  i 


(A-l) 


where 

i(r  =  Temperature  rise  (°C) 

2 

H  =  K/pC  =  Diffusivity  (cm  /s) 

2  0 

K  =  Heat  conductivity  (cal  cm/cm  -5  C) 

2 

p  Density  (gr/cm  ) 

C  =  Heat  capacity  (cal/gr  °C) 

3 

Q  =  Heat-generating  density  (cal/s-cm  ) 
t  =  Time  (s) . 

In  situations  where  heat  is  generated  by  absorption  from  an  electro¬ 
magnetic  wave,  the  heat-generating  density  can  most  conveniently  be  ex- 

2 

pressed  in  terms  of  the  absorption  of  an  incident  fl'ix  P  (watts/cm  ), 
where  the  absorption  by  the  medium  may  be  expressed  in  terms  of  the  ab¬ 
sorption  coefficient. 

If  it  is  assumed  that  a  plane  wave  is  propagating  in  the  z  direction, 
then  the  power  density  in  the  medium  will  be  given  by: 
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dP  o  -ox 
dr  *  1 .  1M 


r<  O 


(cul/s-cm 


(A-2) 


In  n  simpler  model  that  we  will  consider  below,  It  Is  assumed  that 
the  absorbing  layer  absorbs  uniformly.  11'  the  dopth  of  the  layer  Is 
tnken  to  be  2£,  then  the  heat- gene  rating  density  Is  given  by: 


po* 

'1.18  (24) 


f(t) 


(A-3) 


where  0  is  the  fraction  of  the  Incident  power  that  Is  absorbed. 

The  solution  ol  the  heat-conduction  equation  (A-l)  In  three  dimen¬ 
sions  Is  required  to  handle  the  thermal  processes  of  the  retina  during 
exposure.  It  is  assumed  that  the  radiation  is  absorbed  In  small  volumes, 
and  that  the  heat  is  conducted  away  in  all  directions.  Thus,  we  assume 
that  the  heat-generating  volumes  are  located  in  an  infinite  medium  of 
isotropic  properties  and  of  totally  uniform  temperature  as  the  initial 
condition. 

Equation  (A-l)  may  be  rewritten  as: 

=  -uCr“,  t)  (A-4) 


where 


The  Green's  function  for  Eq.  (1)  is 


g(R,T)  = 


R_ 
4ht 


ftySt(rrT) 


3/2 


U(t) 


(A-5) 
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where 


ii  -  \t-  -?1 


t  -j  t  •  i  . 

Thu  solution  of  Eq.  < A- 1 )  is 
-1 


'if?’,  I )  =  ^dl7^  u(i/,l,)GCr‘-“r/,t  - 


,  /  »  ,d  / 

t  )  Cl  !• 


=  I  dr  I  uCr)  l  - 

■'o 


)d  r' 


Using  Eq.  (A— 5 )  in  Eq.  (A-6)  we  obtain,  for  interaction,  volumes 
circularly  symmetric: 


,  2  '2 

,t  r  ®  r 


•liCr^t)  = 


U0  I  dr  4ht  [  1,1 

Jrfc 


4HT  ,  / V 

u  X  f(r  )S  S 


©  z 


0  0 
where  f(r')  is  a  radial  source  function,  and 


S  =  2rr  I 
© 


(— ) 

0\ 2kt  / 


and 


r*  - 

=  I  dz'  e 

J-l 


(z-z ,)2 
4xt 


.  (  z  -  l  z  +  l) 

=  /kttt  [erfc  — p=~  -  erfc  — . _ 

y  |  2/HT  2/h?  I 


Substituting  Eq.  (A-8)  and  Eq.  (A-9)  into  Eq.  (A-7)  we  obtain: 


(A-6) 

that  are 

(A-7) 

(A-8) 

(A-9) 
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# 


(A- 10) 


The  special  case  of  interest  is  the  temperature  along  the  axis--i.e., 
r  =  0.  For  this  case, 


■„(£)  ' 


Wo  consider  now  two  cases: 


Case  I — Uniform  Radial  Distribution. 


f(r') 


1  r  ^  a 


0  r  >  0 


then 


r  ... 


4ht 


=  2kt|1  -  e 


2  1 

t\ 

) 


Substituting  in  Eq.  (A-10)  we  obtain: 

.t 


UK  f 

■■■’  ■  *  / 


Cr*,  t)  =  |  dT  erfc  | -  erfc 
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A 
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C a hu  II—  Gaussian  Itudlal  Distribution. 


f(r') 


then 


1. 


S  /r '  '** 

r'dr'  c  ""  « 


2h~ 


0 

Substituting  into  Eq.  (A-10)  we  obtain: 
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1  *- 


•Ik* 


i(rf,t)  = 
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dr  j  erfc 
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where,  from  Eqs.  (A-3)  and  (A-4), 


( A- 13) 


(A- 14) 
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